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Abstract 

We have investigated some of the properties of dense sub-nuclear matter at the crustal region 
(both the outer crust and the inner crust region) of a magnetar. The relativistic version of Thomas- 
Fermi (TF) model is used in presence of strong quantizing magnetic field for the outer crust matter. 
The compressed matter in the outer crust, which is a crystal of metallic iron, is replaced by a regular 
array of spherically symmetric Wigner-Seitz (WS) cells. In the inner crust region, a mixture of 
iron and heavier neutron rich nuclei along with electrons and free neutrons has been considered. 
Conventional Harrison- Wheeler (HW) and Bethe-Baym-Pethick (BBP) equation of states are used 
for the nuclear mass formula. A lot of significant changes in the characteristic properties of dense 
crustal matter, both at the outer crust and the inner crust, have been observed. 



PACS numbers: 97.60.Jd, 97.60.-s, 75.25.-l-z 



1. INTRODUCTION 



Magnetars are the most exotic stellar objects, believed to be strongly magnetized young 
neutron stars. The surface magnetic field for such objects are observed to be > lO^^G 

innn 

1|, l2|, y, IJ] . Then it is quite possible that the field at the interior, even at the inner crust region 
may be stronger than the surface value (predicted theoretically by scalar Virial theorem). 
If the internal field strength is happened to be so high, then most of the physical and 
chemical properties of dense stellar matter of the magnetars must change significantly from 
the conventional picture . A lot of investigations have already been done on the effect 

of strong quantizing magnetic field on various physical properties of dense stellar matter 
inside neutron stars as well as quark matter inside quark stars or hybrid stars, including 

;he effect on quark-hadron phase transition at the core region of a compact neutron star 

8|. The effect of such strong magnetic field on various elementary processes inside neutron 
stars and quark stars or hybrid stars have also been studied. The effect of strong quantizing 
magnetic field on the /5-equilibration among the constituents have also been investigated 

9j. It has also been shown that strong quantizing magnetic field acts like a catalyst to 
generate fermion mass dynamically, i.e., chiral symmetry breaking occurs in presence of 
strong quantizing magnetic field 10|, lul, Il2 |. 

In this article we have presented our investigation on the effect of strong magnetic field on 
the crustal matter of magnetars. The work is divided into two parts: in the first part, based 

n 

on one of our very recent work p^, we have investigated the effect of strong quantizing 
magnetic field on the outer crust matter and in the second part, we study the properties of 
compact sub-nuclear matter at the inner crust region in presence of such strong quantizing 
magnetic field. 

The paper is organized in the following manner: In section 2, the effect of strong quan- 
tizing magnetic field on the outer crust matter of magnetars is discussed, in section 3, we 
have presented with detailed numerical computation on the effect of strong magnetic field 
on the inner crust matter, which is assumed to be a mixture of iron, some heavier neutron 
rich nuclei, electrons and free neutrons. The presence of free neutrons are considered beyond 
neutron drip density. In this section, we have re-investigated the properties of inner crust 
matter of a typical neutron star in presence of strong quantizing magnetic field. In the inner 
crust region, for metallic iron and more heavier neutron rich nuclei, we have considered the 
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conventional HW and BBP equation of states [is], [3] • Finally, in the last section, we 
have given the conclusions and discussed the importance and future prospects of the present 
work. 

2. OUTER CRUST MATTER 

In a recent work we have developed an exact formalism for the relativistic version of TF 
model in presence of strong quantizing magnetic field [l^. This formalism is used within 
the limitation of TF model [l7| , to obtain the equation of state state for crustal matter of a 
typical magnetar, which is mainly a dense crystal of metallic iron in the sub-nuclear density 
region. In this model, the compressed iron atoms / ions are replaced by the spherically 
symmetric WS cells, with positively charged nuclei at the centre, surrounded by a non- 
uniform cloud of electron gas. In this model it has been assumed that in a WS cell, the 
atomic number which is the number of protons within the nuclei is Z and the mass number 
is A (number of protons and neutrons). To make each cell electrically charge neutral, Z 
must also be the number of electrons inside the WS cells. 

In this model, the modified form of Poisson's equation is given by 

^=X:(2-M(0^(^)-0^a:^)^/^ (1) 

where is the modified form of electrostatic potential related to the original Coulomb 
potential V{x) by the relation 

m = Y^^Ve + ey{x)i (2) 

with /Xe, the electron chemical potential, assumed to be constant throughout the cell (this 
is the so called Thomas- Fermi condition), the dimensionless scaled radial coordinate x is 
related to the actual radial coordinate r of WS cells by the relation: x = r j with 

\ 1/2 



e = |e|, the magnitude of electron charge, B is the constant magnetic field, assumed to be 
along z-direction (we have chosen the gauge = (0, 0, xB, 0)), 

<^o = (4) 
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with = {ml + 2z/ei?)^/^, rrie = O.SMeV, the electron rest mass and z/ = 0, 1, 2, , t'max, 

is the Landau quantum number for the electrons, i^max is the upper limit of Landau quantum 
number (the upper limit of the Landau quantum number u^aax is finite at zero temperature, 
otherwise z/max = oo). Finally, the factor (2 — 6uq) indicates that the zeroth Landau level is 
singly degenerate, whereas, all other states are doubly degenerate. In this article we have 
assumed that the electron gas in the dense crustal matter of metallic iron crystal is strongly 
degenerate and is considered to be at zero temperature. Then it is quite obvious from the 
non-negative value of electron Fermi momentum, that the upper limit of Landau quantum 
number z/max is given by 

To obtain numerical solution for (j){x) for a given magnetic field and a particular set of Z 
and A, we further assume that instead of a point object, the nucleus at the centre of a WS 
cell, has a finite dimension and is assumed to be spherical in nature, so that the corresponding 
radius r„ = tqA^^^, with rg = L12fm. Such a choice also removes the singularity problem of 
TF equation at the origin [l^ . Therefore it is not necessary to follow the prescription given 



by Feynman, Metropolis and Teller to obtain the numerical solution for TF equation 20|. 
Again from the physics point of view, the potential must satisfy the boundary conditions, 
given by 

rV(r) = Ze for r r^, and = for r ^ (6) 

dr 

where is the radius of the WS cell. Then by simple algebraic manipulation it is easy to 
show that the modified form of Coulomb potential satisfies the boundary conditions 

0(a;)U=x„ = l and = (7) 

ax X 

where Xn = rn/ fi, the scaled nuclear radius and Xg = rs/fi, the corresponding scaled radius 
of the WS cell. Both these quantities are dimensionless. 

Further, the right hand side of the Poisson's equation (eqn.(l)) must be real. Which 
requires (poX <\ |. Hence we get an additional condition, to be satisfied by the upper 
limit of Landau quantum number, and is given by 

..„W<(^^W=-^], (8) 

From the definition of Landau quantum number, the above inequality must necessarily be 
> 0. The above equation also shows that the upper limit of Landau quantum number 



depends on the position {x or r coordinates) of the electron within the WS cell, with which 
it is associated. 

Since electron distribution is non-uniform within each cell, the Fermi momentum of a 

particular electron must depend on its positional coordinate in the cell. Then it is expected 
that the variation of pi;'(r) will be such that the electron chemical potential, given by 

/^e = iPF{rf + ml + 2u{r)eBf''^ 

remain constant throughout the cell with proper space dependent Landau quantum number 

V. 

Satisfying all these conditions, of which some of them are particularly necessary for this 
model, we have solved the Poisson's equation numerically within the range of r from nuclear 
surface to the WS cell boundary, for B = Bf^, 10 x S(^), 10^ x Bf^ and 10^ x Bf^, where Bf^ 
is the quantum critical limit for the magnetic field at and above which the Landau levels are 
populated for the relativistic electrons, and is given by eB^^ — with the choice of unit 
Ti = c = \. The magnetic field beyond this limit is called the quantizing magnetic field and 
the quantum mechanical effect plays an important role in this domain. Within the range 

< 2; < Xs, the numerical solution for ^{x) can be fitted exactly by a straight line for a 
particular magnetic field strength B and can be expressed as ^{x) — ax + b. Obviously, the 
parameters a and b are functions of magnetic field strength. In TABLE-I we have shown 
the variation of the parameters a and b with magnetic field strength. 
TABLE-I 





a 


b 


Xs (MeV-i) 


a 




P 


Q 


10° 


-2.48 


1.0002 


30.67 


0.0113 


2.694 


0.04 


0.574 


10^ 


-2.54 


1.0001 


11.78 


0.0789 


2.602 


0.25 


0.591 


10^ 


-2.69 


1.0001 


4.58 


0.7157 


2.504 


0.39 


0.743 


10^ 


-3.08 


1.0002 


1.78 


5.0198 


1.948 


0.52 


0.842 



(Variation of the parameter a and b for — ax + b and scaled surface radius Xg for WS 
cells, parameters a and P for Eke — ckx^ and the parameters p and q for E^^ = px^ with 
magnetic field strength) 

Because of some numerical in-accuracy, the values of the parameter 6, give the intersec- 
tions of the straight lines with y-axis, are not exactly one. Instead of fixing its exact value 
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which is equal to one for all B, and fitting the straight lines with a as unknown parameter, 
we have used our numerical code for general purpose. However, for all the magnetic field 
values, the parameter b is very close to one. The numerical values for the other parameter a 
is always negative and its magnitude increases with the increase in magnetic field B. Which 
actually means that the radius or the volume of a particular WS cell decreases with the 
increase in magnetic field strength. In TABLE-I we have also shown the explicit variation 
of Xs, the scaled surface radius of the WS cells, for various magnetic field strength. While 
solving the Poisson's equation numerically for a given magnetic field strength B, the instruc- 
tion is given in our numerical code to terminate if the surface condition, given by eqn.(7) 
is satisfied and hence wc obtain (also = aj^/i) as a function of magnetic field strength. 
In fig.(l) we have shown the variation of Xg and also the corresponding with the strength 
of magnetic field. This figure clearly shows that the WS cells become more compressed if 
the magnetic field becomes stronger. This is in some sense analogous to what is called the 
magnetostriction in classical magneto-statics. These two curves can also be fitted by the 
power law functions, given by 



Knowing the scaled Coulomb potential (f){x) at various x points (in the range x„ to Xg) within 
the WS cell for a given magnetic field strength, we have evaluated numerically Vmaxix) at 
every x points within the cell. The variations of i^max(2;) with x for four different magnetic 
field strengths arc shown in fig. (2). Although Vnn^xix) must be a set of discrete numbers, for 
the sake of illustration, we have plotted it as a continuous variable. In this figure, curves a 
and b are for B — B^^ and 10 x B^") respectively. For these two curves, the variables J^max(a^) 
is plotted along the left side y-axis and x (this is actually x^ in Mev~^) along x-axis at the 
bottom. Similarly for S = 50 x Bf^ and 100 x B^f> , the variations are shown by the curves 
a and /3 respectively. In this case, i^max(a^) is plotted along the right side ?/-axis and x is 
plotted along upper x-axis. From this figure, it is possible to make a number of conchisions: 
(i) For larger B values, Xg are smaller, (ii) For smaller 5, Pui3x{x) starts with quite large 
value near the nuclear surface, e.g., = 124 and = 15 as shown in curves a and b respectively. 
Whereas, for curve a it starts with i/max = 2 and for /3, the starting value is i^max = 1- (hi) 
The discrete nature of v-aiax is obvious from the high magnetic field curves a and (5. (iv) 
Finally, for all the field values, the upper limit fraaxi^) becomes exactly zero at the surface 




MeV 



-1 



and 
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of the WS cells. In other wards, we can say that all the electrons near the WS cell surface 
are strongly polarized and the spins are anti-parallel to the direction of magnetic field. This 
is, of course, a purely relativistic effect. The possibility of fully polarized scenario can easily 
be obtained from the analytical solution of Dirac equation for electrons in presence of strong 
quantizing magnetic field. The eigen functions will not be simple spinor solutions, whereas 
the energy eigen value will be = [pi + 2veB + m'lY^'^^ with 2i/ = n + 1 + m^, where 

= 0, 1,2, .. is the Landau principal quantum number and mg = ±1, the eigen values for 
the spin operator Q, 21]. Hence, for u = Uja&x = 0, the only possible choice is the 
combination n = and iris = — 1- Which actually means that in the zeroth Landau level the 
spins of all the electrons are in the direction opposite (this is due to negative charge carried 
by the electrons) to the external magnetic field. We have further noticed that beyond the 
field value 100 x B^'> (which is > lO^^G), the upper limit t'max(a;) becomes identically zero 
not only at the surface region of WS cells, but at all the points inside the cell. 

To obtain the density distribution of electrons within the WS cells, let us consider the 
expression for electron number density, given by 

ne{x) = ^^ Y: {2-6^o)Pf{x) (10) 

where the electron Fermi momentum Pf{x) can be obtained from the TF condition and is 
given by 



Mx) = <!^V(^^j -ml} (11) 

Since the electron density is larger near the nuclear surface, i.e., at the central region of the 
WS cells compared to the values near their boundary regions, the effect of electron density 
dominates over the influence of magnetic field strength at the central region, whereas near 
the WS surface, since the density is low enough, it will be of completely opposite picture; the 
magnetic field will play the significant role. This will make the upper limit i^max large enough 
(except for B > lO^^G) near the nuclear surface and quite small or identically zero near WS 
cell boundary. The results are also quite obvious because the electron Fermi momentum is 
minimum at the nuclear surface and maximum near the the WS cell boundary. Later we shall 
see that the electron kinetic energy part will also behave like Fermi momentum. Therefore, 
the quantum mechanical effect of low and moderate values for quantizing magnetic field will 
be more significant near the surface region of WS cells. Near the core region, the effect is 
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almost classical. This is particularly valid for the curves a and h of fig. (2), because for these 
curves z/max(2;) started with quite large values. 

It is also obvious from eqn(ll), that the electron Fermi momentum is a function of the 
positional coordinate of the associated electron. Using the numerically fitted form of 
we obtain Pf{x) for a given magnetic field strength and the corresponding ne{x) for the 
same magnetic field. In fig. (3) we have plotted the variation of electron density within WS 
cells for four different magnetic field strengths. In this figure, the curves indicated by the 
numbers 0, 1, 2 and 3 are for B = B^^\ = 10 x B'f \ B = W x B'f^ and S = 10^ x B'f^ 
respectively. It is quite obvious from these curves that the electron density increases with 
the increase in magnetic field strength. Further, for all the values of 5, electron density is 
maximum at the centre and minimum at the surface. 

Let us now evaluate the variation of (i) electron kinetic energy, (ii) electron-nucleus 
interaction potential, (iii) electron-electron direct interaction potential and (iv) electron- 
electron exchange interaction within the WS cells. The kinetic energy part for electrons at 
a particular point (r) within the WS cell is given by 

Eke{x)= / (2-M / dpM + <f"-^e] (12) 

Evaluating the integral over analytically and then substituting for Pf{x) from eqn.(ll), 
with the fitted form of (f){x), one can obtain the electron kinetic energy Eke, as a function 
of r or X by the numerical evaluation of above integral over r or x respectively. We have 
seen that within the WS cells, for a particular value of B, the kinetic energy part satisfies a 
power law of the form E^e = cnx^, where the parameters a and P are functions of magnetic 
field strength. In TABLE-I we have shown these variations. The change is more significant 
for a than (3. Therefore, just like the electron density, the electron kinetic energy is also an 
increasing function of radial coordinate x. 



Next we consider the three possible types o 
the electron- nucleus interaction part, given by 



interaction potential. Let us first consider 
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r 



-A-nZe^jj^ / xdxn(.(x) (13) 

J Xri. 



To obtain the electron- nucleus interaction energy within the WS cell for a constant 5, 
we substitute the expression for electron number density ne(a;) from eqn.(lO), and then 



numerically integrate over x In fig. (4) we have shown the variation of Een{x) = \Een{x)\ 
as a function of x within the WS cell for three different magnetic field strengths: B = 
B'f\ W X Bf> and 10^ x B'f^ For very low x values, i.e., almost on the surface of the 
nucleus at the centre, the magnitude of the potential energy is an increasing function of 
X, which actually means that the region is strongly attractive in nature. Again at the 
surface region, just at the skin of the cell, it is again an increasing function of x (with larger 
gradient), so that the attractive field again becomes extremely strong at the surface region 
to keep the outer most electrons confined within the cell. At the middle region, there is 
a kind of saturation and force due to electron-nucleus interaction almost vanishes. This is 
quite analogous to the normal metallic scenario, in which electric field can not exist. 

Next we consider the electron-electron interaction. Let us first evaluate the direct term. 
It is given by [l^ 

= ie^l dhn,{r) J rfVne(r )p^ (14) 

Assuming r as the principal axis and 6 is the angle between r and r', we have d^r = Anr^dr, 
d^r' = 27Tr'^dr' sin 6d6 (we have assumed that the vectors r and r' are on the same plane) 
and |r — r'l = (r^ + r'^ — 2rr'cos6')^/^. The limits for both r and r' are from r„ to and 
the range of 6 is from to vr. Let us first evaluate the angular integral, given by 

sin 6 



I{r,r')= r 
Jo 



(r^ _l_ j-'2 _ 2^^i ^)i/2 

It is straight forward to show that I{r,r') = (r + r') — \r — r'\. Then we can express the 
electron-electron direct part in the form 

E^f = 47r^e^ T rdrne{r) f r' dr'ne{r')[{r + r') - \r - r'\] (15) 

Now from this equation it is trivial to show that the quantity within the third bracket will 
be 2r' for r' < r and for the opposite case it will be 2r. Then the above expression for direct 
interaction, reduces to the following simple form: 



E^^) = Se^TT^I / rdrn^{r) / r''^dr'n^{r') 

+ / r'^drne{r) / r'dr'ne{r')^ (16) 

Following the same procedure adopted for e — n case, we have evaluated numerically the 
above coupled integrals. In fig. (5) we have shown the variation of E^^{x) within the WS cell 



for three different magnetic field strengths: B = 10^ x B^^ and 10^ x B^\ Surprisingly, 
the variations are almost identical with \Een\- However, the magnitudes are several orders 
less than the corresponding electron-nucleus part. Of course, unlike the E^^i the direct term 
is positive throughout and for all the values of magnetic field strength. 

Next we shall consider the exchange term, which is negative in nature. The magnitude 
of exchange energy integral corresponding to the iih. electron in the cell is given by [3] 

^Ir^ = ^J2 [ c^VrfV— i^V^,(f)V^,(fO^,(r-')^,(fO (17) 

J I' ' I 

where the spinor wave function iplf) is given by eqns.(2)-(5) in [l^ and ^/^(r) = '?/'"''(r)7o, the 
adjoint of the spinor and 70 is the zeroth part of the Dirac gamma matrices 7^. Now it is 
very easy to show that for t = t' 

— —* '2tti 

Mr)ipi{r') = exp[-i{py{y - y') + p,{z - z')}] 

)} (18) 

Similarly, we have 

iljj{r')ijj{f) = exp[i{p'y{y - y') + p'^{z - z')}] 

{Iu';pdx)Iu';pdx') + 1^' ^l.p'U) L' ~l;pdx')} (19) 



where Iu;py{x) is same as ly, given by eqn.(5) in [l^. When these two terms are combined, 
we have, after replacing the sum over j by the integrals 

n+OO r+PF 



/+00 r+PF 
dp'y / dp'^ 
-oo ■ J—PF 



ee 



2 / V LILIE,, 



^(2-5,,o) / ... / LydpJyLJpJ^d\d\'—-^^ 
ex^[-i{{py - p'y){y - y') + [p, - p'^){z - z')}] 

[{Iu;py{x)lu;py{x') + I^-l^p^x) Iu-l;Py{x')} 

Uu';pdx)Iy>,pdx') + Iu'-l;pdx)Iu'-l;pdx')}] (20) 



It is possible to evaluate the integrals over y' and z', given by 22 1 



— oo J — oo 



hoo r+oo ]_ 

47r 

= _exp(-ir|x-x'|) (21) 
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where K — [{py — py^Y + [p^ — Pz'YY^'^- Further, the integral over y and z is given by 



Since the final form of exchange integral is mult i- dimensional in nature with an extremely 
complicated structure of integrand, it is absolutely impossible to have analytical solution for 
e — e exchange. We have therefore adopted the multi-dimensional Monte-Carlo integration 
code to evaluate the exchange term. It is found the variation of the magnitude of exchange 
energy with the scaled radius x within the WS cells for a constant B can be expressed by 
the power law formula, given by 



where the parameters p and q are again functions of magnetic field strength B. The variation 
of the parameters p and q with magnetic field are shown in TABLE-I. 

From the variation of the parameters p and q with 5, it is quite clear that the magnitude 
of exchange energy also increases with x and also with the magnetic field strength B, i.e., 
minimum at the central nuclear surface and maximum at the WS cell boundary. 

The total energy of an electron at any point x within the cell is then given by (with 
negative values for exchange energy) 



We have seen that at the central part up to certain radial distance Tc from the centre of the 
WS cells, the total energy is negative, whereas beyond this radial distance and up to the 
WS cell surface, it is positive. We have evaluated numerically the x {— Xc say), the value of 
the point within the cell at which the total energy of the electrons just become zero for a 

given magnetic field. The position Xc of this zero energy point changes with the strength of 
magnetic field and is shown in fig. (6). Further, this variation can also be expressed by the 
functional form, given by 



with a = 0.5, (5 = 1.34 and 7 = 0.205. If the total energy for the electrons at the central 
part is found to be negative, then obviously, we can not expect that all the Z electrons in 
the cell particularly near the central region, are participating in statistical process. Some of 




Eif=px\ 



Etot = Eke + Een + E^f + E^f . 




(22) 
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them remain bound in negative energy states near the nuclear surface at the centre of the 
WS cell. 

The overall pressure term can also be obtained for the electron gas from the total energy 

Etot, given by 



where 



P{x) = nl{xf-^ = P,(x) + PXx) (23) 
one. 



P,._g(2-M/^ (24) 



is the electron kinetic pressure part. The momentum integral for the kinetic pressure can 
very easily be obtained, and is given by 

_|_ 



P.(p^ + m^)^/--m-lnf ^^ + ^^^ + --^^ 



(25) 



with pf — pf{x). This expression gives the electron kinetic pressure at various points within 
the WS cell. The other term, P^., coming from the interaction part, Ec — -Bgn + + E^^^\ 
and is given by 

BE 

Pc{x) = nlix)^^, (26) 

The interaction part of electron pressure is obtained numerically by substituting the ex- 
pression for Pf{x) and evaluating the derivatives over Ue numerically for a given S, at a 
particular point x within the cell and also for a given set of A and Z. We have noticed that 
at the central region of WS cells, since the density dominates over the magnetic field part, 
the total pressure is positive, while at the outer end, the surface region, since magnetic field 
dominates over the density contribution, it becomes negative. It is found numerically that 
at some point within the WS cell for a given B, the total pressure becomes exactly zero. 
It is found that this critical position is a function of magnetic field strength and decreases 
with the increase in magnetic field strength. In fig. (7) we have plotted this variation. The 
functional dependence can also be expressed as 

B 



a + P exp 



-7 



MeV-^ (27) 



with a — 0.38, P — 1.28 and 7 = 0.171. The qualitative form of this variation is exactly 
identical with that of Xc for total energy. 
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3. INNER CRUST MATTER 



In this section we shall study the effect of strong quantizing magnetic field on dense 
nuclear matter composed of nuclei, free neutrons and electrons at sub-nuclear density. In 
the conventional neutron star model, it is generally believed that for a typical neutron star 
of radius ~ 10km, the width of inner crust is about 0.5 — 0.6km. The constituents are a 
mixture of heavy nuclei (of which some of them are highly neutron rich), electron gas and 
there are also free neutrons (if the matter density is above the neutron drip point). It is 
also believed that the free neutrons in the inner crust region may be in super-fluid state. 
The density range of this zone depends on the type of equation of state or the mass formula 
for the heavy nuclei, including the highly neutron rich nuclei. The density range in the 
conventional picture is ~ 10^ — lO^^gm/cm^. The inner crust is covered by the outer crust 
and is mainly the crystalline structure of compressed solid crystal of fully ionized metallic 
iron immersed in a dense electron gas. In the previous section we have made a detailed 
study of various aspects of this sub-nuclear matter in the outer crust region. For a typical 
neutron star, the width of the outer crust region is ~ 0.2 — 0.4Km. Analogous to the model 
calculation of the outer crust matter, which is assumed to be a regularly arranged structure 
of WS cells, for the study of inner crust region, we have assumed that the matter is a mixture 
of WS cells of not only metallic iron, but a large number of neutron rich nuclei and a gas of 
free neutrons. We further assume that the free neutrons in this region, above the neutron 
drip point are in normal fluid state. The effect of strong quantizing magnetic field on the 
super-fluidity of free neutron matter, in which magnetic field interacts with the anomalous 
magnetic moment of neutrons, will be studied in future. 

In this section, to study the effect of strong quantizing magnetic field on inner crust 
matter, we consider two types of conventional nuclear mass formulae, which are generally 
used below the nuclear saturation density, but applicable near neutron drip region. These 
two mass formulae are (i) Harrison- Wheeler (HW) equation of state, and (ii) Bethe-Baym- 
Pethick (BBP) equation of state Q, Q, 16|. 



In HW equation of state, the relativistic electrons are assumed to be in /3-equilibrium 
with the nuclei (which also includes the neutron rich nuclei). Whereas, above the neutron 
drip density, the /5-equilibrium is among the nuclei, free neutrons and the electrons. In this 
article, we assume /5-equilibrium configuration in presence of a strong quantizing back ground 
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magnetic field. In this section, the properties of inner crust matter have been investigated 
in presence of such a strong magnetic field. In the conventional physical picture, in presence 
of free electron gas in the medium, the balance between the Coulomb force and the nuclear 
force, which gives Fe*^ as the most stable nucleus will get shifted towards the heavier nuclei. 
The nuclei, formed in such an environment will contain more neutrons (by inverse /3-decay) 
compared to the usual picture. The Coulomb force plays a very little role. The nuclei 
becomes more and more neutron rich as the electron density goes up and when the matter 
density increases to ~ 4 x lO^^gm/cm^ in the inner crust region, the ratio n/p reaches a 
critical level. Any further increase in density will lead to neutron drip in the medium. In 
such a scenario, highly neutron rich nuclei, electrons and free neutrons co-exist in chemical 
equilibrium. Again, with the increase in matter density, beyond neutron drip density, free 
neutrons appear in the medium and at some stage the kinetic pressure of the system will 
be dominated by the free neutron pressure. As we shall see that this positive pressure 
contribution of free neutron gas will make the total pressure within the inner crust region 
a positive definite, even in presence of strong magnetic field. Whereas in the outer crust, 
the total pressure, coming from the electron gas only and is negative at the surface region 
of WS cells in presence of strong magnetic field. Just like the outer crust region, we have 
assumed that only electron part of inner crust gets affected by the strong magnetic field. We 
have not considered the effect of quantizing magnetic field on bound protons and neutrons 
and above neutron drip density, on free neutrons. We have noticed that the presence of 
strong magnetic field makes a lot of qualitative and quantitative changes in the physical and 
chemical properties of inner crust matter. 

The BBP equation of state is applicable for nuclear matter in the density range from 
neutron drip density pdrip to normal nuclear density pnuc- Here also the matter is assumed 
to be composed of nuclei [Fe^^ and other neutron rich nuclei beyond Fe^^), electrons and 
free neutrons. It is well known that BBP equation of state is a considerable improvement 
over the HW equation of state. In BBP equation of state, a mass formula, more or less 
like the HW equation of state is used but incorporated a lot of improvements from detailed 
many body calculation. The nuclear surface energy, for example, assumed in the previous 
treatment to be that of a nucleus in vacuum. An introduction of free neutron gas out side the 
nuclei reduces the nuclear surface energy. This is quite correct, because when inside and out 
side of a nuclei become identical, the surface energy must vanish. In BBP equation of state. 
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the nuclear Coulomb energy is included in a more accurate manner and is called nuclear 
lattice Coulomb energy. It is well known that the BBP equation of state is applicable up 
to nuclear density, therefore, when all the bound nuclei dissolves into a continuous matter, 
mainly composed of free neutrons and a tiny fraction of protons and electrons, to incorporate 
this type of melting process of nuclei, in the BBP equation of state the factors giving the 
fractional volume occupied by the nuclei and the fractional volume occupied by the free 
neutron gas have been taken into account. It is also assumed that at this density (< pnuc.) 
the nuclei are stable to /3-decay. Further, the neutrons in the free neutron gas are in chemical 
equilibrium with the electrons and the nucleons within the nuclei and in addition to the f3- 
stability of the nuclei, the whole system must be in chemical equilibrium. In this model the 
kinetic pressure of free neutrons must necessarily be equal with the pressure of the nucleons 
bound within the nuclei. This is the condition for mechanical equilibrium. In our present 
calculation, with this BBP equation of state, we have investigated some of the interesting 
properties of neutron star inner crust matter in presence of strong quantizing magnetic field, 
assuming that only the electron part gets affected by the presence of strong magnetic field. 
Since detailed mathematical formalism for both HW and BBP equation of states, which 
are not affected by strong magnetic fields, are available in a large number of classic papers, 
including the original ones, and also included in many standard text books jl^ . [isl . 16|. 



we are therefore do not feel the need to include those mathematical steps in this article. 
Further, the effect of strong quantizing magnetic field on the electron part has already been 
discussed in the previous section. In this section, we therefore present only the numerical 
results based on these two equation of states and the results related to electron gas as given 
in the previous section. 

In the numerical computation, we have found that for very low mass number, since there 
is no free neutrons available in the system, even much above the neutron drip density, only 
electrons contribute in total kinetic pressure and is negative. This result is almost identical 
with the outer crust scenario. In fig. (8) we have plotted the variation of the critical value of 
mass number A and the corresponding atomic number Z with the strength of magnetic field, 
using HW equation of state for the nuclei and for the electrons, for which the resukts are 
given in the previous section. These are the critical values at which the kinetic pressure of the 
inner crust region just becomes zero, i.e., the system just becomes mechanically stable. It is 
obvious from this figure that the minimum of mass number of the nuclei for which the inner 
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crust matter just becomes mechanically stable, increases with the strength of magnetic field. 
In other words, this is equivalent to say that the presence of strong magnetic field makes 
the nuclei more and more massive or increases the minimum mass of the nuclei in the inner 
crust region. In some of our previous work, done long ago we have seen that the same 
conclusion is applicable for quark matter. The presence of strong quantizing magnetic field 
generates quark mass dynamically in dense quark matter composed of massless quarks. This 
is the well known magnetic field induced chiral symmetry violation. The strong quantizing 
magnetic field acts like a catalyst to generate / increase mass. In fig. (9) we have plotted the 
same kind of variation using the BBP equation of state. In this case the minimum mass of the 
nuclei are more than that obtained from HW equation of state. Further, one should notice 
from these two figures (figs. (8) and (9)) that for low and moderate strength of magnetic field, 
the effect is not so significant. But beyond lO^^G, when all the electrons occupy the zeroth 
Landau level, or in other words, when the quantum mechanical effect of strong magnetic field 
is most important, the minimum mass rises sharply to very large values. One can see from 
12] that the qualitative nature of the curve showing the dependence of dynamical quark 
mass on the magnetic field strength is exactly identical, although the physical scenarios are 
completely different. 

In fig. (10), we have plotted the variation of the ratio Ue/n and n„/n with the mass 
number using HW equation of state. Here rig is the electron density, is the free neutron 
density and n = Un + n^AjZ is the total baryon density of inner crust matter. Dashed 
curve is for 5 = 0, middle and the lower curves, indicated by el and e2 are for B = lO^S^^^ 
and lO^i?^^-* respectively. The curve indicated by n is for free neutron gas. In the HW 
equation of state, the neutron number density above the neutron drip point is independent 
of magnetic field strength. Further, the mass number around which neutrons are liberated 
from neutron rich nuclei is about 95, which is again independent of magnetic field strength. 
However, immediately after the emission of neutrons from heavy neutron rich nuclei the 
overall kinetic pressure can not become non-zero. This is also quite clear from the figs. (8) 
and (9). In fig. (11), we have plotted the same kind of variations as shown in fig. (10). Here 
we have used the BBP equation of state. Solid curves are for B = 0, while the dashed curves 
are for m^B^^^\ In this case the neutron drip out from heavy neutron rich nuclei around 
A = 100, which is little bit heavier than the HW case. Further, the qualitative nature of 
Xn = nn/n is totally different from HW case. Instead of increase initially and then saturates. 
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as we observe in HW case, it decreases monotonically with A. Also, the free neutron density 
very weakly depends on the strength of magnetic field. These qualitative differences are 
found to be the consequence of chemical equilibrium among the free neutrons, electrons and 
nucleons within the heavy neutron rich nuclei and the overall /3-equilibrium condition. 

In fig. (12). the variation of total nuclear density is plotted against the mass number 
from HW equation of state for four different magnetic field strengths:i? — 0, B — 1QB^'^\ 
B = lO^S^^) and B = 10^B^^\ indicated in the diagram by 0, 1, 10^ and 10^ respectively. 
Initially, for low A values, only the bound neuclons within the nuclei (= Arie/Z) contributes. 
Just beyond the critical value, A — 95, since neutrons drip out from the heavy neutron rich 
nuclei, the number density jumps suddenly. Of course, below A = 60, we can not have stable 
inner crust matter. In fig. (13), we have plotted the same kind of variation for BBP equation 
of state. In this case the qualitative and the quantitative changes are again because of the 
chemical equilibrium conditions 

The equation of state from HW mass formula is plotted in fig. (14). Here the lower curve 
is for B — 0, middle and upper curves are for B — Bf^ and B — l^fB^^ respectively. In 
fig. (15), we have plotted the same kind variation from BBP equation of state. 

From these two figures it is quite obvious that in presence of strong quantizing magnetic 
field, the inner crust matter becomes mechanically stable only at very high density, when 
enough number of free neutrons are available in the system. In this high density situation, 
positive neutron kinetic pressure, dominates over negative electron pressure, in presence of 
strong quantizing magnetic field. This is true for both HW and BBP equation of states. 
Further, it is also obvious from the figures, that the matter becomes softer with the increase 
in magnetic field strength. The qualitative nature of both the equation of states arc almost 
identical in presence of strong magnetic field. However, the HW equation of state is a bit 
softer than the BBP case. 

4. CONCLUSIONS 

In the first part of this article we have presented our investigation on the properties of 

dense outer crust matter of the magnetars. We have replaced the dense metallic iron crystal 
by a regular array of spherically symmetric WS cells. It has been observed that the radius 
of each cell decreases with the magnetic field strength. In this article, however, we have 
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not considered the magnetic field induced deformation of WS cells, which is quite hkely if 
the magnetic field is strong enough. In future we shall present this important issue with a 
completely different approach. 

It has also been noticed that the upper limit of Landau quantum number is a function 
of positional coordinate of the electron with which it is associated within the WS cells. We 
have observed that at the surface region, for all the values of magnetic field strength, this 
upper limit becomes identically zero. Which actually means that the electrons near the 
WS cell surface are strongly polarized in the opposite direction of external magnetic field. 
Whereas, for B > lO^^G, they are polarized at every points within the cells. 

It has been observed that the electron density within the cells increases with the increase 
in magnetic field strength. Further, for all the values of magnetic field strength, the electron 
density is maximum near the nuclear surface (r = r„) and minimum at the WS cell boundary 
(r = Vs). 

We have also studied the variations of kinetic energy, electron-nucleus interaction energy, 
electron-electron direct potential energy and electron-electron exchange interaction part. We 
have shown the variation of these quantities within the WS cells for a number of magnetic 
field strengths. Interestingly, we have noticed that the total energy of the electrons is 
negative near the nuclear surface within the cells and becomes positive beyond some radial 
distance (xc) which increases with the strength of magnetic field at low and moderate region 
but saturates to some constant value for B > lO^^G, when all the electrons occupy their 
zeroth Landau level. Hence we can not conclude that all the Z electrons in the cell are 
participating in statistical process. Some of them remain bound in negative energy orbitals 
near the nuclear surface at the centre of the WS cell. 

In the second part of this work, we have investigated some of the properties of inner crust 
matter of magnetars. We have used HW and BBP equation of states for the nuclear mass 
formula. It has been observed from both the models, that for a stable inner crust matter, 
the nuclei present must be heavier than iron and much more neutron rich. The heaviness is 
more in the case of BBP equation of state. We have also noticed that for low and moderate 
values of magnetic field strength, the variation of mass number and the corresponding atomic 
number with magnetic field is not so significant. Whereas, for B > lO^^G, when electrons 
occupy only the zeroth Landau level, then much more heavier neutron rich nuclei are formed 
in the inner crust region. It is found that high magnetic field behaves like a catalyst to 



18 



generates heavy neutron rich nuclei. 

We have observed that initially the electron density increases with the increase in mass 
number, but as soon as free neutrons appear in the system, the electron density decreases 
and saturates to some constant value which depends on the magnetic field strength. In 
the case of HW equation of state, free neutron density does not depend on the strength of 
magnetic field, whereas, for BBP case, because of chemical equilibrium condition, the free 
neutron density depends very weakly on the magnetic field strength. We have noticed that 
in the case of BBP equation of state the overall qualitative difference is because of chemical 
equilibrium among the constituents. 

The total baryon density rises sharply like an avalanche for the value of A at which free 
neutrons appear in the system. However, for BBP equation of state, because of chemical 
equilibrium condition, the rise is not so sharply visible for a given magnetic field B. 

The qualitative nature of equation of states are almost identical. It is found that in 
presence of strong magnetic field, the inner crust matter becomes mechanically stable (with 
the positive value of kinetic pressure) only at very high density. 

We therefore believe that to study various properties of dense matter associated with 
magnetars, one must consider all these significant changes from the conventional scenario. 
Further, we expect that the strong magnetic field, if present well within the inner crust of 
magnetars, must affect the super-fluidity of cold neutron matter. 
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FIG. 1: The variation of scaled surface radius Xg in MeV ^ (solid curve) and the actual radius 
in A(dashed curve) with the strength of magnetic field strength (expressed in terms of B^^K 
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FIG. 2: The variation of Umax with x. For the curves a and b, the relevant axes are at the bottom 
(x) and at the left side (z^max), whereas for the curves a and /?, these axes are at the top and at 
the right side respectively. 
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FIG. 3: The variation of electron density in MeV within the cells for B = h x Be with h 
1, 10, 100 and 1000 indicated by the numbers 0, 1, 2, and 3. 
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FIG. 4: The variation of the magnitude of electron-nucleus interaction energy within the WS cells 
ioT B = hx B'f^ with /i = 1, 100 and 1000 
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FIG. 5: The variation of electron-electron direct interaction energy within the WS cells for B = 
h X B^""^ with /i = 1, 100 and 1000 
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FIG. 6: The variation of critical position (xc) within the cell at which Etot is exactly zero with the 
strength of magnetic field (* is the numerically evaluated points and f is the fitted point). 
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FIG. 7: The variation of critical position (xc) within the cell at which total pressure is exactly zero 
with the strength of magnetic field (* is the numerically evaluated points and f is the fitted point). 
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FIG. 8: Variation of critical values for A and Z with the strength of magnetic field B (expressed as 
B/B'f^) at which kinetic pressure becomes just zero. This result based on HW equation of state 



200 




FIG. 9: Variation of critical values for A and Z with the strength of magnetic field B (expressed as 
B/ B^f^) at which kinetic pressure becomes just zero. This result based on BBP equation of state 
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FIG. 10: Variation of the ratio ng/n and n„/n with the mass number using HW equation of state. 
Here Ue is the electron density, n„ is the free neutron density and n = n„ + ricA/Z is the total 
baryon density of inner crust matter. Dashed curve is for B = 0, middle and the lower curves are 
for B = 10^ Sc^-* and lO^Bc^^ indicated by el and e2 respectively. The curve indicated by n is for 
rin/n. 
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FIG. 11: Variation of the ratio Xg = ng/n and x„ = n„/n with the mass number A using BBP 
equation of state. Here Ue is the electron density, n„ is the free neutron density and n = Un+rieA/Z 
is the total baryon density of inner crust matter. Solid curve is for B = and the dashed curve is 
for 10 X Be Electron part and the neutron part are indicated by Xe and Xn respectively. 
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FIG. 12: Variation of total baryon density with the mass number of the nuclei, using HW equation 
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of state for four different magnetic field strengths:^ = 0, B = 105c , B = IQ Be and B = 
10^ 5^""^ indicated by 0, 1, 10^ and 10^ respectively. 
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FIG. 13: Same kind of variation as in fig. 12 but using BBP equation of state. 
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FIG. 14: Equation of state from HW mass formula. Lowe curve is for B = 0, middle and upper 
curves are for B = B'f^ and B = lO^B^^^ respectively. 
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FIG. 15: Same kind of plot as shown in fig.l4 but for BBP equation of state 
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